Cooper pair tunneling in circuits with substantial dissipation: 
the three-junction R-pump for single Cooper pairs 
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We propose a circuit (we call it R-pump) comprising a linear array of three small-capacitance 
superconducting tunnel junctions with miniature resistors (R > Rq = h/Ae 2 = 6.5 kQ) attached to 
the ends of this array. Owing to the Coulomb blockade effect and the effect of dissipative environment 
on the supercurrent, this circuit enables the gate-controlled transfer of individual Cooper pairs. The 
first experiment on operating the R-pump is described. 
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INTRODUCTION 

The Coulomb blockade effect in circuits with small- 
capacitance tunnel junctions provides the means of ma- 
nipulating single charge quanta (see, for example, the 
review paper by Averin and Likharev jl)). If periodic 
signals of frequency / are applied to the gate electrodes 
of the circuit, a train of single charges q can move across 
an array of junctions so that charge pumping, giving rise 
to the current / = qf, is achieved. It was experimentally 
proven that the normal-state metallic circuits enable sin- 
gle electrons (q = e) to be effectively pumped at frequen- 
cies / of about several MHz ||. Moreover, the accuracy 
of single-electron pumping can nowadays meet the re- 
quirements of fundamental metrology, viz. SI / 1 ~ 10~ 8 

Unlike pumping of electrons, the pumping of Cooper 
pairs (g = 2e) in superconducting circuits has not been 
that successful so far. The only experiment had been 
carried out in 1991 by Geerligs et al. with a three- 
junction Al sample. Although this experiment did evi- 
dence a pumping of the pairs, the pumping was strongly 
disturbed by several factors: the Landau-Zener transi- 
tions, Cooper pair co-tunneling, quasiparticle tunneling, 
etc. As a result, the shape of the current plateau at 
I = 2ef in the I-V curve was far from being perfect. 

Recently, in their theoretical paper Pekola et al. || 
concluded that pumping of Cooper pairs with reason- 
able accuracy in a three-junction array was impossi- 
ble. For practical values of parameters (the ratio of 
the Josephson coupling energy to the charging energy 
A = Ej/E c = 0.01 — 0.1) they evaluated the inaccu- 
racy of pumping 61 /I to be as much as 9%-63%. This is 
because of intensive co-tunneling of pairs in the short ar- 
rays, i.e. the process of tunneling simultaneously across 
several (> 2) junctions. To suppress the co-tunneling 
and, in doing so, to improve the characteristics of the 
pump they proposed to considerably increase the num- 
ber of junctions 7V(^> 3) and, hence, of the gates (N— 1). 
However, such modification would make operation of the 
circuit more complex. 

In this paper we propose an alternative way to im- 
prove Cooper pair pumping in a three-junction array. 
We modify the bare circuit by attaching to the ends of 



the array the miniature resistors (see the electric dia- 
gram of this device, which we call R-pump, in Fig. |l|). 
Their total resistance R exceeds the resistance quantum 
Rq = h/Ae 2 = 6.5 k£! so that the dimensionless parame- 
ter 



R 
Rq 
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The self-capacitance of these resistors should not be much 
larger than the junction capacitance C . In this paper we 
will analyze how such resistors affect tunneling and co- 
tunncling of Cooper pairs in the three-junction array and 
present preliminary experimental data. 



PECULIARITIES IN OPERATING THE R-PUMP 

In general terms, the principle of operation of the R- 
pump remains the same as that of the pump without 
resistors (for details of the Cooper pair pump operation 
see Ref. Q). Two periodic signals V\{t) and ^(t) are 
applied to the gates to form an elliptic trajectory in the 
parameter plane V\ — Vi. At zero voltage V across the 
pump, this trajectory encircles in the clockwise (counter- 
Islands (Al) 
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FIG. 1: Schematic of the R-pump for single Cooper pairs. 
Three superconducting junctions of type Al/AlOx/Al form 
two small islands. Miniature normal metal (Cr) resistors ac- 
complish the structure. The device is driven by two harmonic 
voltages Vi,2 = Foi, 02 + ^4cos(27r/i ± i#), where the dc offset 
Foi, 02 determines the centre of the cycling trajectory in the 
Vi — Vi plane; the phase shift 6 > 
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clockwise) direction the triple point of the boundaries 
between the stability domains for adjacent charge config- 
urations with an excess Cooper pair. A complete cycle 
results in the charge q = 2e (q — — 2e) being sequentially 
transferred across all three junctions of the array. 

The Cooper pair transitions in a circuit without dis- 
sipation (quasiparticlc tunneling is neglected) occur due 
to Josephson tunneling (supercurrent). The strength of 
the Josephson coupling Ej = fi-I c (here I c is the critical 
current) is assumed to be -C E c = e 2 /2C, the charg- 
ing energy, i.e. parameter A C 1. When the domain 
boundary is crossed, the resonance condition is obeyed 
and Josephson coupling mixes up the two states, each 
corresponding to the position of the pair on either side 
of the junction. (Note that, at the same point, the reso- 
nance condition for co-tunneling in the opposite direction 
across two other junctions is met.) Slow passage across 
the boundary results in transition of the pair in the de- 
sired direction. This process is mapped on the adiabatic 
level-crossing dynamics (see, for example, Ref. Q), and 
the probability of missing the transition is given by the 
Landau-Zener expression: 



Plz = exp 
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where E oc f is the velocity of variation of the energy 
difference E between the two states. 

In contrast to the elastic (and, hence, reversible) tun- 
neling of pairs in the pump without resistors, tunneling in 
the R-pump is accompanied by dissipation of energy. The 
strength of the dissipation effect depends on the active 
part KeZ(uj) of an electromagnetic impedance seen by 
the tunneling pair, namely by the dimensionlcss parame- 
ter z' = ReZ(O) /Rq. If the condition of weak Josephson 
coupling Ej -C E c /\fzJ is satisfied, i.e. \\f~z' <C 1, the 
rate of Cooper pair tunneling is expressed as |?j 



(3) 



where E is the energy gain associated with this transition 
and P(E) is the function describing the property of the 
environment to absorb the energy released in a tunneling 
event. (The discussion of properties as well as numerous 
examples of P{E) for different types of electromagnetic 
environment are given in the review paper by Ingold and 
Nazarov ||.) The probability of a missed transition is 
then given by the expression: 



exp 



(4) 



Here we used the normalization condition 
J_ P{E)dE — 1, and the assumption of a value 

E constant in time was exploited. Thus, in the case 



of weak coupling, dissipation does not affect the total 
probability of tunneling p = 1 — p m . This conclusion is 
in accordance with the general result obtained by Ao 
and Rammer |jj for quantum dynamics of a two-level 
system in the presence of substantial dissipation. 

Using the expansion of P in the region of small E -C E c 
0, H one arrives at the expression 



T(E) ex X 2 E 
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For z' < A the rate T is peaked at E — 0, while in the 
case (z' > the rate is zero at E = and the maximum 
is reached at finite energy E < 4E Cl or, in other terms, at 
finite positive voltage V across the junction. This prop- 
erty of T(E) in circuits with substantial dissipation plays 
the crucial role in the operation of the three-junction R- 
pump. 

For one junction in the three-junction array with gate 
capacitances C g <§C C the effective impedance ReZ(O) = 
R' = i.e. the parameter z' — The factor 

i, which considerably attenuates the effect of resistors, 
stems from the square of the ratio of the total capacitance 
of the series array (= -|C) to the junction capacitance C 
(see a similar network analysis in Ref. || ) . 

Accordingly, the equivalent damping for two junctions 
of the three-junction array is larger and determined by 
R' = |i? and, hence, z' — |z. Therefore, the rate of 
Cooper pair transition simultaneously across two junc- 
tions, i.e. the rate of the co-tunneling, at small E is 



T cot {E) oc X 4 E 
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This expression has been obtained by considering the en- 
vironmental effect in a similar fashion as was done by 
Odintsov et al. |J for the co-tunneling of normal elec- 
trons in a single-electron transistor. Accordingly, co- 
tunneling across all three junctions decays even stronger 
(oc X 6 E 2z ~ 5 ), because the whole array experiences the 
full resistance R. 

The comparison of Eq. (6) with Eq. (5) taken for the 
value z' = ^z shows that the rate of co-tunneling is dras- 
tically depressed if z is sufficiently large. For example, 
for z — 9 (or R w 58.5 kf2) the rate of Cooper pair 
co-tunneling across two junctions r cot tx E 5 (i.e. it is 
similarly small as single-electron co-tunneling across a 
three-junction normal array Q). On the other hand, the 
direct tunneling rate T oc E, is similar to that of ordi- 
nary electron tunneling. The effective junction resistance 
in the latter case is R e s ~ 0.03 A -2 Rq. For the practi- 
cal value A ~ 0.03 this formula yields R c g ~ 200 kfi, 
i.e. the typical value of a normal single electron junc- 
tion. Since the junction resistance determines the char- 
acteristic time constant RC , it can be expected that the 
frequency characteristics of the Cooper pair R-pump are 
at least not worse than those of the normal electron coun- 
terpart without resistors. 
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EXPERIMENT 

The Al/AlOz/Al tunnel junctions and Cr resistors 
were fabricated in situ by the three-angle shadow evap- 
oration technique through a trilayer mask patterned by 
e-bcam lithography and reactive-ion etching. The tunnel 
junction parameters of the best suited available sample 
were found to be: capacitance C m 250 aF and normal- 
state resistance Rj w 160 kf2. The former value yields 
E c w 320 /ieV. On the assumption that the Ambegaokar- 
Baratoff relation between the critical current I c and the 
junction resistance is valid, the aforementioned value of 
Rj yields the value Ej = Aa\Rq /2Rj ~ 4.1 ueV, where 
Aai ~ 200 £ieV is the superconducting energy gap of 
aluminum at low temperature. These parameters of the 
sample give the ratio of characteristic energies A « 0.013. 

The thin-film Cr resistors (thickness d w 7 nm) had 
lateral dimensions of about 80 nm by 10 fim. For the 
purpose of an independent characterization of the resis- 
tors and junctions, we attached four resistors in twos to 
either end of the junction array, resulting in the equiv- 
alent resistance \R on each side as depicted in Fig. [j], 
(For the measurement of the I-V curves we used, how- 
ever, the two- wire configuration.) Either of four resistors 
had R w 60 kfl and reasonably low self-capacitance per 
unit length, w 60 aF//im. The resistors exhibited negli- 
gibly small (< 2 %) non-linearity of the I-V curve in the 
region of small voltages at the millikelvin temperature 




The samples were measured in the dilution refrigera- 
tor at the temperature T of about 10 mK, i.e. well be- 
low the critical temperature of aluminum, T c w 1.15 K. 
The bias and gate lines were equipped with pieces of the 
thermocoax cable which was thermally anchored at the 
sample holder plate and served as a filter for frequencies 
above 1 GHz Q. The typical I-V curves with and with- 
out ac drive are presented in Fig. |2[ The curves exhibit 
considerably smeared steps whose position (I — 2ef) 
clearly shows the linear current-versus-frequency depen- 
dence (although the quality of the steps is noticeably de- 
graded with increasing frequency) fll2| . Another feature 
of the steps is their position on the voltage axes: they ap- 
pear at finite voltage V applied in the current direction, 
while at V = neither co-tunneling nor pumping effect 
were observed. We attribute this effect to the damp- 
ing effect of the resistors, which was too strong in this 
sample. Thus, the behavior of the Cooper pair R-pump 
differs significantly from that of the Cooper pair pump 
without resistors Q. 

DISCUSSION 

Although our first experiment evidenced the desirable 
effect of Cooper pair pumping, a further improvement is 
inevitably needed. First, the sample parameters should 
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FIG. 2: I-V characteristics of the R-pump without and with ac 
drive of two different frequencies. The inset shows the blow- 
up of the steps whose positions were expected to be strictly 
at / = — 2ef (shown by dashed lines), where the sign "-" is 
due to the phase relation chosen between two gate signals. 



be optimized: Josephson coupling should be increased 
(possibly up to A = 0.1) while the resistance value R 
should be somewhat reduced. These parameter values 
lead to an increase of the tunneling rate and, as a conse- 
quence, to Cooper pair pumping at zero voltage bias V. 
Secondly, the quasiparticle tunneling, resulting in spo- 
radic translations of the cycle trajectory and, hence, to 
the pumping errors, should be reliably depressed: As the 
measurements showed, all our Al-Cr samples (as well as 
those with E c < Aai) so far suffered from quasiparticle 
tunneling JD|] and did not show the so-called parity effect 
fhij , |ll| . It is particularly remarkable that our Cr resis- 
tors do not prevent an uncontrolled poisoning of Joseph- 
son tunneling by non-equilibrium quasiparticles arriving 
at the islands from the external circuit. Such "buffer" 
effect of the normal electrodes was proposed and demon- 
strated in Al-Cu devices by Joyez et al. |l(|. That is why 
the problem of the poisoning quasiparticle tunneling re- 
quires a more radical solution. Probably manufacture 
of small-capacitance niobium junctions (ANb ~ 1-4 meV 
3> E c ) could improve the situation by making the parity 
effect in these structures stable. 

In conclusion, we proposed a simple superconductor- 
normal metal circuit enabling, in principle, the efficient 
pumping of Cooper pairs. Due to the energy dissipation 
in the resistors, the effect of Cooper pair co-tunneling 
is heavily damped. A possibility of pair pumping was 
demonstrated and improvements of the experiment are 
proposed. 

This work is supported in part by the EU Project 
COUNT. 
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